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1.  IVIRODUCTION 

Observed  wind  and  temperature  profiles  used 
as  basic  states  in  models  often  contain 
modifications  by  the  secondary  flow,  the  correct 
basic  state  however,  is  the  one  existing  prior  to 
the  onset  of  roll  development.  A  low-order 
spectral  model  allows  for  efficient  study  of 
secondary  instability  and  its  effects  on  the 
initial  basic  state  variables.  Using  an  Eknan 
wind  profile  as  an  idealised  case,  we  show  that 
errors  in  the  Initial  wind  profile  can  lead  to 
significant  errors  in  the  model  results, 
particularly  in  vadues  of  the  preferred  roll 
orientation  angle.  Ihe  ultimate  application  of 
the  otodel  however,  is  to  observed  botndary  layers. 

2.  MODEL  SUttARY 

A  nonlinear  14-coefficient  spectral  model 
of  thermally  and  dynamically  forced  shallow 
Boussinesg  convection  is  derived  using  truncated 
Fourier  expansions  for  the  stream  function  t  and 
perturbation  temperature  T  (Haack  and  Shirer 
1991).  Iheae  expansions  are  ccxaposed  of  two 
parts:  roll  circulation  patterns  given  by  Tr  and 
^ri  snd  background  modification  profiles  given  by 
Tb  and  d%. 

The  roll  patterns  are  specified  by  four 
components  in  the  Fourier  expansions  that  involve 
a  horizontal  wavenuaber  of  one  and  vertical 
wavenumbers  qsl  and  t%=2  (Haack  and  Shirer  1991). 
Wo  also  specify  four  components  In  the 
perturbation  tesqperature  expansion  that  represent 
the  background  temperature  modification: 

Tb=T|  sin(  qr  z/Zpj.)  4T]ain  ( nrs/z^ ) 

♦TjsinC  (q-n)re/2;j,]«T4sin(  (q+n)»z/Zj,I 

(1) 

and  two  in  the  stream  function  expansion  for  the 
backgtxxsid  wind  modification: 

l^(z)sU|Cos{  (q-n)rz/Zp] 

•fUioos{(q+n)rz/z^]  (2) 

The  model  accepts  as  input  an  observed  or 
idealized  wind  profile,  U)  and  V, ,  in  standard 
east/north  coordinates  that  is  rotated  into  a 
cross-roll  wind  U(z)  by  introducing  a  roll 
orientation  angle  fii 

U(z)sU«sin4  -  V,cos4  (3) 

Baoauss  the  oirculationa  are  asstned  to  be 
two-disHnsional  in  struoturs  and  the  (k>riolis 


tenaa  are  neglected,  the  along-roll  wind  V(z)  is 
eliminated  from  the  system.  Thus,  the  basic  state 
wind  profile  is  represented  in  the  model  by  the 
following  expression: 

U(  z )  =Uo '•■Ui  ain(  qrz/Zpj,)  ■t^U200s(  qrz/z^) 

Hl3sin(nrz/2;^)-fU4COs(nrz/z^)  (4) 

involving  five  Fourier  coefficients,  two  vertical 
wavenuBibers  q  atnd  n,  and  the  domain  height  z.^,.  In 

general,  we  find  that  (4)  gives  an  excellent 
representation  of  most  atmospheric  wind  profiles 
when  q=l  and  n=2.  In  addition  to  the  Fourier 
coefficients,  the  dynamic  forcing  is  also  given  in 
the  dimensionless  system  by  the  Reynolds  nimber 


(6) 

where  |V(z^)|  is  the  wind  speed  at  the  domain  top 
and  )f  is  the  constant  eddy  viscosity. 

The  basic  state  tesqwrature  is  more  simply 
specified  by  a  linear  profile  obtained  from  the 
vertical  gradient  AsT.  This  quantity  includes  the 
combined  effect  of  the  environment^  lapse  rate 
and  a  vertically  distributed  surface  buoyaiwy 
contribution.  The  corresponding  thermodynamic 
forcing  in  the  dimensionless  system  in  defined  by 

the  Rayleigh  nimber 

gz^sT 

Ras— 2- -  (6) 

where  c  is  the  constant  thenDometrio  conductivity 
and  Toe  is  the  surface  air  teoiperature. 

The  model  geometry  is  illustrated  in  Fig.  1 
for  a  standard  Eknan  wind  profile  having  a 
westerly  geostrophic  wind  (Haack  and  Shirer 

1991).  The  y-axis  is  aligned  parallel  to  the  rolls 
that  are  propagating  slowly  northwaud  in  the  -x 
direction.  A  positive  roll  orientation  angle  fi  is 
measured  oointerclockwise  from  the  standard 
direction  x«  (east)  to  the  roll  axis  y.  The 
critical  values  Ree  and  Rsc  of  the  forcing  rates, 
idiich  are  foizid  by  a  stability  analysis  of  the 
conductive  solution,  depend  on  both  the 
orientation  angle  t  and  the  aspect  ratio  a:2z^L, 

where  L  is  the  horizontal  wavelength.  The  values 
of  Rse  and  Rsc,  when  minimized  with  respect  to 
both  /  and  a  (or  L) ,  give  the  preferred  or 


s 


figure  1,  Sekemeiic  depietioe  of  »  fropegoliog  roll  streem  fuMciioB 

fttere  gioeo  kg  He  14-coeffieie»t  model  of  Meeck  e»d  Skirer  (1991). 
Mere  *  smell  posilive  orieeletioa  eagle  fip  is  skoma  for  e  case  milk  e 
ateaderd  Skmea  profile  tket  kes  e  mesterly  geostropkic  miad  vector 

Ike  exes  (xt,Pt)  represeal  Ike  ateaderd  (eeat/aorlk)  coordiaete 
spatem,  ead  Ike  exes  (x,p)  repreaeat  tke  roll  coordiaete  apstem,  Ike 
boaaderp  leper  deptk  is  deaoled  bp  tj  ead  ike  roll  apeciag  bp  Lp. 

toll  propegetioa  is  domamiad  from  rigkl  to  left. 


expected  values  of  the  roll  tfeoaetry.  These 
expected  values  are  labeled  4p  and  Up  in  Fig.  i 
and  are  the  ones  coopared  with  observations. 

3.  RESULTS 

Shown  in  Fig.  2  is  the  curve  of  critical 
forcing  values  for  the  transition  fron  a 

conductive  solution  to  a  periodic  roll  solution  as 
given  by  the  linear  stability  analysis.  Here  we 
use  an  Eknan  depth  of  D=190  ■,  which  gives  model 
results  for  Re<60  that  agree  with  those  produced 
by  a  higher  resolution  nusarical  analysis  of  the 
linear  model  (Haack  and  Shirer  1991).  Labeled 
along  the  curve  are  the  preferred  values  of  4 
a.  The  path  arrows  specify  a  possible  transition 
to  the  supercritical  forcing  rates  indicated  at 
points  (A)  and  (B) .  At  these  values  of  the 

forcing  rates,  we  integrate  the  model  and  examine 
the  roll  patterns  and  modifications  to  the  basic 
state. 


The  integrations  are  carried  out  until 
stable,  energetically  steady  roll  solutions  are 
obtained.  This  generally  requires  less  than  four 
hours  so  that  the  neglect  of  the  Coriolis  terms  is 
a  valid  assunption.  For  integration  point  (B) , 
the  roll  structure  in  dimensionless  quantities  T* 

and  d*  are  shown  in  Figs.  3a  and  3b.  For  points 

(A)  emd  (B),  profiles  of  the  background 
modification  in  dimensional  quantities  Tb  and  Ubi 
given  by  (1)  and  (2),  are  shown  in  Figs.  4a  and 
4b.  These  profiles  suggest  that  the  roll 
modifications  act  to  reduce  the  thermal  and 
dynamic  instabilities  that  initially  produced  the 
secondary  flow.  The  basic  temperature  is  edtered 
by  cooliiig  the  lower  layer,  autd  warning  the  upper 
layer  which  produces  a  near-neutral 
stratification.  However,  the  magnitude  of  this 
modification  is  not  likely  to  be  measurable  in  the 
atmosphere  nor  great  enough  to  significantly 
affect  the  stability  results  shown  in  Fig.  2. 

However,  as  shown  in  Fig.  4b,  the  roll 
modifications  to  the  basic  wind  profile  can  be 
significant.  Here,  the  rolls  alter  the  background 
wind  shear  by  as  much  as  1.2  ms*  1/200  n  in  the 
center  of  the  domain.  Using  observed  winds  as 


figure  t.  Treaaitioa  curve  repreaeatiag  tke 

miaimum  values  of  critical  foreiag  aeeded 
for  roll  developmeat  mkea  usiag  e  ateaderd 
Skmea  miad  profile  keviag  Skmea  deptk 
9*1 90m.  fref erred  values  of  roll 
orieatelioa  eagle  4p  ead  aspect  ratio  Op 
ere  labeled  eloag  tke  curve.  Tke  petk 
erroms  ladtcele  ea  etmospkeric  treasilioa 
to  tke  supercritical  foreiag  rates  at 
poiats  i  ead  t. 


input  to  the  model,  we  have  found  modifications  of 
up  to  a  few  meters  per  second  (Shirer  and  Kaack 
1990).  Fis.  5  shows  the  original  and  the  altered 
cross-roll  Ekman  profiles.  The  rolls  create  a 
dynamically  more  stable  environment  by  reducing 
the  cross-roll  shear. 

In  Fig.  6,  a  comparison  of  the  original 
Ekman  spiral  and  the  modified  one  for  integration 
point  (A)  illustrates  the  characteristics  of  the 
altered  flow.  We  note  that  the  turning  angle  in 
the  modified  Ekman  spiral  is  reduced,  and  the 
altered  flow  has  maximun  velocities  and  is  more 
closely  aligned  with  the  direction  of  the 
geostrophic  wind  at  low  levels.  These  results  are 
consistent  with  those  found  by  Brown  (1970)  whose 
model  of  secondary  circulations  forced  by  Scman 
flow  had  Coriolis  forcing  terms. 

To  examine  the  consequences  of  using 
roll-modified  winds  in  boundary  layer  models,  we 
use  the  point  (B)  modified  cross-roll  wind  profile 
shown  in  Fig.  5  in  the  stability  analysis.  The 
resulting  transition  curves  and  preferred  xcll 
(dtaracteristics  a  and  4  for  the  modified  wind  are 
shown  in  Fig.  7.  Upon  oomparing  these  results 
with  those  in  Fig.  2,  we  find  that  the  transition 
to  a  roll  solution  occurs  only  in  the  statically 
unstable  regime  (Ra>0)  because  of  the  reduced  wind 
shear  in  the  modified  profile.  Moreover,  values 
of  preferred  roll  characteristics  in  Fig.  2  differ 
markedly  from  those  in  Fig  7.  Significantly, 
values  of  4p  are  altered  by  a  maxiimmi  of  14*  at 
low  values  of  Re.  This  error  in  /)p  is  of  the 
order  reported  by  Shirer  and  Brumaer  (1986)  in 
which  the  contributions  by  the  rolls  were  not 
considered. 

4.  OONCUBIONS 

A  14-coefficient  spectral  model  of  shallow 
Boussinesq  flow  was  used  to  study  modifications  of 
initial  bsckgrouixl  wind  and  temperature  profiles 
by  boundary  layer  roll  circulations.  Temporal 
integrations  at  supercritical  values  of  the 
dynamic  and  thermal  forcing  rates  Re  and  Ra 
produced  stable  periodic  roll  solutions  and 
time-independent  nonlinear  wind  and  temperature 
modification  profiles.  Alterations  by  the  rolls 
of  the  basic  temperature  were  dimensionally  small, 
while  those  of  the  basic  wind  created  modified 
profiles  with  significantly  reduced  cross-roll 
shear.  When  a  modified  Ekman  wind  spiral  was  used 
in  a  linear  stability  analysis,  roll  orientation 
angles  shifted  by  a  maximum  of  14*  from  those 
given  by  an  analysis  using  the  original  winds. 
Oxisequently,  roll  perturbations  can  markedly 
change  the  structure  of  backgrowd  wind  profiles 
—  even  in  the  absence  of  Coriolis  forcing.  This 
secondary  flow  contribution  should  be  removed  from 
wind  measurements  before  using  the  observed 
profiles  as  basic  states  in  boundary  layer  gxxlels. 

5.  ACXNOWLEIXMENTS 

This  research  has  been  partially  supported 
by  the  National  Science  Foundation  through  Grants 
ATM86-198S4  and  A‘m89-09947,  the  National 
Aeronautics  and  Space  Administration  through 
Contract  NAS8-361S0  and  Grant  NAC8-780,  and  the 
Office  of  Naval  Research  through  Contract 
N0000t4-86-K-06880.  NOARL  contribution  no. 
91:111:432.  Approved  for  public  release; 
distribution  is  unlimited. 


figure  S.  Peiierus  of  dimeusiouleas  roll 

temferature  ia  (a)  aai  atream  fuaetioa 

d*  is  (h)  for  foiai  B  ia  fig.  i. 


figure  4-  frofilea  of  iimeaaioaal  roll 

temperature  modificatioa  ia  (a)  aai  roll 

miad  modificatioa  u^  ia  (b)  for  poiata  d 

aai  i  ia  fig.  t. 


U  +  u^  (m/s) 


figurt  5.  Ike  origieel  cross-rell  Skmee  mini 
profile  (solid)  lokee  ,  ted  the 
roll-modified  Ekmtt  profiles  for  poitis  i 
(dotted)  ted  i  (dtsked)  it  fig. 2. 
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